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OPTIMUM FLIGHTS WITH PRESSURE OF LIGHT 
TAKEN INTO ACCOUNT 


by 

S. N. Klrpichnikov 


This article has for its subject the construction of single-impulse tra- 
jectories of interorbital flights with a minimum consumption of mass in the 
gravitational field of the Sun — a spherical central body radiating light symmet- 
rically. It disregards the perturbing action on the space vehicle caused by 
objects in space located on the specified boundaries of the Kepler orbits. It 
assumes that all three orbits, the initial, intermediate, and final, are located in 
one plane and the motion along these orbits is proceeding in one direction. 

The pressure of light can be, on one hand, regarded as a small correc- 
tion factor for flights of conventional spaceships. On the other hand, the pres- 
sure of light becomes comparable with the force of the Sun’s gravity in cases of 
research flights of probes, which are made of thin hollow shells (cylinders) 
whose inside is filled with low-pressure gas and whose outside is covered with 
a good light- reflecting coating. The use of such shells was suggested by Ehricke 1 
for investigating the cosmic space of our solar system and also for transporting 
useful loads. It notes the simplicity of these probes, their light weight, and 
their ability to carry useful loads. Such objects have an intense brightness, 
which makes it possible to increase by many times the probability of successful 
observation of the probes with the aid of telescopes. These shell-probes can be 
launched from the Earth and also directly from spaceships or from artificial 
satellites. 


1 K. Ehricke, "Instrumented Comets-Astronautics of Solar and Planetary 
Probes, " VIII International Astronautical Congress, Barcelona, 1957. 
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1. The Pressure of Light 


Let us find the main vector and the main moment of forces of light pres- 
sure acting upon a stationary body whose radiated surface is S. The derivation 
of expressions in this chapter will be based on the quantum theory of light; the 
effect due to anisotropism of the reradiation will be neglected. 

Let us introduce a stationary system of coordinates x, y, z; the z-axis 
coincides with the heliocentric radius-vector of the irradiated body; i. e. , it is 
directed along the parallel beam falling on the body. Let us separate an elemen- 
tary area (dS) on the body's surface, with R as its coefficient of reflection. 

Let p = {£, 77, t,} be the radius-vector of some point on (dS) , the origin of this 
vector is located at point 0 of the irradiated body; a is the angle of incidence 
equal to the angle of reflection; p is the angle between the projection of the out- 
side normal to the area on plane xy and the x-axis counted in positive direction 
from the x-axis. In this case, the elementary change in momentum during the 
time dt, will be 


dK = (M 2 cei - Mjcejjdt , (1) 

where c is the velocity of light; §1 is a unit or basis vector in direction of the 
incident beam; e 2 is a unit-vector in direction of the reflected beam of rays; 

Mj, M 2 are, respectively, the mass of the photons incident to the area (dS) and 
reflected from it per unit of time. 


Based on the principle of equivalence of mass and energy, let us find 
E E 

Mj = — cos a dS, M 2 = R — cos a dS . (2) 

c 2 c 2 

Here, E is the solar constant (intensity of solar radiation per unit of area) for 
the area located at a distance r from the Sun, calculated as: 


E = 



E^r+ = 0.302 • 10 33 ergs/sec , 


(3) 


where r+ is the average distance between the Earth and the Sun and E+ is the solar 
o o 

constant for the distance r+. Using the theorem of momentum, we find the force 
of the light-pressure f on the area ( dS) 
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E \ 

f = -R - sin 2x cos a cos (3 dS , I 
X c I 


f = -R — sin 2 a cos a sin 0 dS , 

y c 


f = — (1 + R cos 2a) cos a dS . I 
z c / 


(4) 


It is then easy to find the projection of the main vector F and of the main 
moment L for the point 0 on the axes of the coordinates 


F = 


F = 


F = 


- — f f R cos 2 a sin a cos (3 dS , ^ 

C (S) 


- — f f R cos 2 a sin a sin 0 dS , 
° (S) 


(5) 


- J f (1 + R cos 2a) cos a dS , / 

0 (S) 


L = — J'Jfrjfl + R cos 2a) + t,R sin 2 a sin /?] cos a dS , \ 
x c <S) 


L = - — f f [£,R sin 2a cos /3 + ^(1 + R cos 2a) ]cos a dS , 
y C (S) 


\ ( 6 ) 


L = — f f R(tj cos p - | sin p) sin 2a cos a dS . 
2 c <S) 


Let us find the pressure exerted by light on a spherical body with radius 
R and with an index of reflection R identical for the entire surface. Let us 
introduce on the body's surface a spherical system of coordinates (longi- 
tude <p, latitude i>) with the origin in the center of the sphere; the longitude cp 
will be counted from the direction to the Sun; therefore, 

cos a = cos <p cos ip . (7) 

Due to the symmetry, the system of forces of light pressure is reduced 
to a resultant force applied to the geometrical center of the sphere and equal to 
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ffR 2 . 


(8) 


J- - o 

z cr z 

Consequently, the force acting upon a stationary, spherical, irradiated 
body is directed along its heliocentric radius-vector and is independent of the 
reflecting capacity of the latter; this result is fully similar to the result of 
Radziyevskiy. 1 Therefore, it makes sense when a well-reflecting coating is 
used only to prevent an excessive overheating. 

NOTE 1. If the irradiated body is moving with a velocity U, both the 
value of the force of light pressure and its direction will change by a value of 
about U/c (for more details, see Radziyevskiy*) . In view of the insignificance 
of the last ratio, this change will be disregarded. 


2. Heliocentric Trajectories with Pressure of Light 
Taken into Consideration 


Let us consider a body with a mass m moving in fields of Newton gravita- 
tion and of the Sun's radiated light, whose mass will be designated by M. The 
components and of the forces of the pressure of light will be disregarded. 

These components are exactly equal to zero for spherical bodies. 


The force of attraction to the Sun is k 2 Mm/r 2 , where k 2 is the gravitational 
constant and the force of repulsion by light is B/r 2 which, for a sphere with 
radius R, is 


B = 


V* 


ttR 2 


(9) 


while B for other bodies can be readily found by using the expressions ( 3) and 
( 5) . Composing the equations for the relative motion of the irradiated body and 
assuming that m « M, we find for its heliocentric radius-vector r 


« k 2 M _ B _ 
r = - — — r + — - r . 

r 3 mr 


(10) 


Let us introduce the "reduced” mass of the Sun 


1 V. V. Radziyevskiy, "Braking by Radiation in the Solar System and the 

Growth of Rings of Saturn," Astronomicheskiy Zhurnal (Astronomical Journal) , 

29, 3, 1952. 
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M' = M( 1 - 5) , 


( 11 ) 


where <5 is a parameter characterizing the "decrease" in the mass of the Sun and 
is equal to 


5 = 


B 

k 2 mM 


( 12 ) 


Consequently, with the "reduced" mass of the Sun, all expressions in a problem 
for two bodies will be valid and, formally, it is only necessary to replace 
everywhere the value 


K = WI 


(13) 


with a new value of 


K = K\/l - 5 (14) 

We will assume that the parameter is <5 elO, 1) . At 6 = 1, the gravita- 
tional attraction will be balanced by the repulsion of light and rectilinear inertial 
flights are possible in any direction. At <5 > 1, the irradiated bodies will fly out 
the solar system along hyperbolas in which the Sun is located in the external 
focus. 


Let us estimate the value of <5 for thin hollow spherical shell-probes hav- 
ing a radius R, a shell -thickness h, and a density y . Assuming that the shell- 
probe carries a useful load whose mass is m 0 , therefore, the overall mass of 
the probe, is 


m = m 0 + 47rR 2 yh 


(15) 


and the parameter <5 will be 


v !» r2 Hi 

K 2 c(m 0 + 4irR 2 yh) K 2 C ^ + 

If m 0 = 0, then 8 is independent of the radius of the probe and is dependent only 
on the product yh. Let us calculate the values of 8 for different h and at 
y = 1 gram/cm 3 , by assuming that m 0 is much smaller than the mass of the 
shell: 
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6 0.5 0.1 0.05 0.01 

h, microns 0.38 1.9 3.8 19.0 


0.005 0.001 

38.0 190 


This estimate is qualitatively in agreement with the estimate of Ehricke 


l 


The constant elements of the orbits with "reduced" mass of the Sun will 
be called geometrical. We will use for such elements: 


P = 


7T 


q ~ 4T ’ T ’ 


(17) 


where I is a focal parameter, e is the eccentricity, co is the angular distance of 
the pericenter, and T is the moment of passing through the pericenter. For the 
independent variable which determines the position in the orbit we will use the 
polar angle i* , so that the positive direction of the counting coincides with the 
direction of the motion. 


The values pertaining to an osculating orbit will be provided with a sub- 
script "ock." Therefore, we have: 


U r = Kq ock Sln (* " W ock) = KN/1 ' 6q sin (^ - o>) , (18) 

’ K [ P ock + q ock 008 (* ' "ock)] " K ‘ JTrT IP + q cos ( » - «> ], (19) 


1 2 
— = rr 


P ock + P ock q ock C0S (* - W ock) = P + pq cos ’ < 20 > 


i*2 - CO 


ock 


du 


ip = K(t 2 - t j) = / . N . 

n J p , (p , + q , cos uV 
t?i - co . ock ( ock ock / 

1 ock x 


! <!/ 


f 

J d(d + 


du 


n/1 - 6 „ J p(p + q cos u) 2 ’ 


( 21 ) 


where tj , t 2 are, respectively, the moments of time corresponding to t*i , t> 2 ; U , 
U are, respectively, the radial and transversal velocity components. Note 


1 K. Ehricke, op. cit . 
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that in the first integral it is necessary to take into account only the explicit 
relationship of the integrand function of u. Both integrals in (21) are readily 
taken ( for example , for an elliptical motion through an eccentric anomaly: the 
osculating for the first integral and the geometrical for the second) . From the 
relationships (18) to (21) we can readily find 


P ock n/T - 5 ’ 


%ck 



- 2pqS cos (t? - to) 


1 - 6 ) 


( 22 ) 


(23) 


tgco 


_ q( 1 - 6) sin co - p6 sin £ 
ock q( 1 - 6) cos w - p6 cos $ 


i? - co 


T , = T + 
ock 


du 


1 v w 

r 

Kn/1 - 5 ' p(p + q cos u) 2 


i? - co 


ock 


du 


- f 

K „ p ock( p ock + Vk cos u )‘ 


(24) 


(25) 


here, the numerator and denominator signs in expression (24) coincide with the 
signs of sin go ^ and cos respectively; in expression (25) , the osculating 

elements should be considered to be constant for calculating the second integral. 
Note that all osculating elements ( during an elliptical motion in geometrical 
elements) are periodic functions of the angle J with a period 2ir. 


3. Setting Up the Problem of an Optimum (in Terms 
of Energy) Flight with the Effect 
of Pressure of Light Taken into Account. 

The System of Necessary Conditions. 

It is necessary to accomplish a flight of a space vehicle within given 
boundary orbits with the aid of a single thrust — an impulse at the initial orbit. 
The value of the characteristic velocity of this impulse is reduced to a minimum 
to assure a minimum consumption of fuel. All orbits are coplanary. The pres- 
sure of light is taken into account only for the intermediate orbit. Depending on 
the specific physical problem, the space vehicle may be the shell-probe described 
above, or any other spaceship for which the pressure of light must be taken into 
account. 
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Let the p ls q t , oji , and Tj be the elements of the initial orbit and p 2 , q 2 , 
co 2 , and T 2 be the elements of the final orbit. It is assumed that the ratio of the 
effective transverse section to the mass of the objects in space is small at these 
orbits and the elements of the orbits are practically unaffected by the pressure 
of light. 

During the moment of time t 1( when the polar angle is equal to d-\, there 
takes place the thrust-impulse which causes the space vehicle to change its flight 
to the intermediate orbit. During the moment of time t 2 when the polar angle is 
the intermediate orbit intersects the final orbit and the space vehicle collides 
with an object in space located on the final orbit. 

With the elements of the boundary orbits specified, it is necessary to 
find the intermediate orbit of the flight, i. e. , to find its geometrical elements 
p, q, and co; also the angles ^ 2 and the moments of time t t , t 2 which will 
reduce to a minimum the characteristic velocities of the initial impulse. After 
these unknown values are found, we will have for the moment of passage through 
the per ic enter 


i?. - co 


1 ir 

T = t i-i57rrr / 


du 

p(p + q cos u) 


7 , i = 1, 2 


(26) 


and the expressions (22) to (25) can be used to calculate the osculating elements. 

The following conditions must be satisfied: 

<P\ = P 2 + pq cos (i?i - co) - pi - pjqt cos (i?t - u>i) = 0 , (27) 

<Pi = P 2 + PQ cos (i? 2 - to) - pi - p 2 q 2 cos U 2 - to 2 ) = 0 , (28) 

<p 3 = ip t + 4> + ip 2 - a = 0 , (29) 

where the function of ip is expressed through the geometrical elements as per the 
expression (21) , and 


"CO. 

7 7 


du 


.(p. + q. cos u) 2 ’ 


i = 1, 2 , (30) 


a = K(T 2 - TO 


(31) 
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The relationships (27) and (28) indicate a continuity of the radii-vectors at the 
starting and finishing points, while the relationship (29) represents the condi- 
tion for the coincidence of the time spent in motion prior to reaching the final 
point, on one hand, along the initial and intermediate orbits and, on the other 
hand, along the final orbit. 

The characteristic velocity AU of the initial thrust can be expressed 
through the elements of the initial and intermediate orbits 1 for similar computa- 
tions) as follows: 


AU = K {q 2 ( 1-6) - p 2 ( 1 - 6) + q? + 3p? - 2p?6 - 2p ^ N/1 ~ ~ 


-2qq 1 N/l - 6 cos (coi - co) 


- 2qjN/l - 6 ^p 


P? 


- Pi 


2-6 1 

1 - 5 J 




COS 


(^ ~ ^l)l 


V* 


(32) 


and for the angle of inclination of thrust $ ( which is counted from the transversal 
in a direction opposite to the motion) , we have 


tg$ 


q^ 1 - 1 sin (i?i - 0 )) - qi sin (<fr - coi) 

( — — l - — l) [ Pl + qi cos (^ - COi) ] 

P 


(33) 


where the signs of the numerator and denominator coincide, respectively, with 
the signs of sin $ and cos $. 


Let us find the minimum of the function 

(AU) 2 

S ~ 2KW1 - 6 


(34) 


in a class of variables p, q, go, t*i, t *2 which are dependent and related with the 
conditions ( 27) to (29) , i. e. , we find ourselves in a class of a conditional 
extremum of a function of a finite number of variables. Introducing the constant 
factors \ 2 , X 3 , we compose a Lagrange function 


1 S. N. Kirpichnikov, "Optimum Complanary Flight Between Orbits," 
Vestnik LGU ( Leningrad University Herald) , No. 1, 1964. 
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g + Z Vi • 

i - 1 


(35) 


As we know, the partial derivatives of Lagrange functions must be equal to zero 
for all variables 


Jjr + t • ^t] sin ( * - “■> 
inU ' -“ )1+ - f^) = 

MPj'fe sin U, - - pq sin (J, - o>) 1 + = 


+ Ai [piqi sin (^ - cui) - pq sin 


(36) 


-p\fi - 6 + 


o , (37) 

cos (Ji - Wj) + \j [2p + q cos (th - u>) ] 

(38) 


dip 

+ X 2 [2p + q COS (t?2 - w) ] + X 3 “ = 0 , 

dp 


, dip 

q\/l _ 5 _ qj cos (u>i - u>) + Xjp cos (t?i - co) + X 2 p cos (^2 " w) + X 3j ^- = ® » 

(39) 


dip 

-qqi sin (u>i - co) + Xipq sin (t?i - co) + X 2 pq sin (i? 2 - u>) + X 3 -g^j = 0 • 

(40) 

The obtained equations (36) to (40) together with the equations (27) to 
(29) form a system of required conditions made up of eight equations containing 
eight unknown p, q, w, t?i, 1 ^, Xj, X 2 , X 3 . These equations must be solved 
jointly. 


NOTE 2. If a problem is under consideration in which specific motions 
are not taken into account, i. e. , the initial configuration of objects in space is 
arbitrary, the condition (29) for the coincidence of the time spent in motion 
prior to reaching the final point should be omitted, while in the remaining system 
of seven equations it is necessary to assume that X 3 = 0. 

NOTE 3. The pressure of light can be readily taken into account at the 
final orbit. For this is sufficient to consider the elements p 2 , q 2 , a T 2 as 
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V 

geometrical and to replace in all equations the function of ipz with = ip 2 ( 1 - 52) 2 

where 6 2 is the parameter 6 computed for an object in space located at the final 
orbit. 


4. Flight Between Circular Orbits 

Assume that the initial and final orbits are circular with radii and r 2 , 
respectively, and their values are qi = 0 and q 2 = 0. 

It can be shown that from the system of required conditions will follow 

that 


\ 3 = 0 . (41) 

Therefore, it is first necessary to solve the problem without taking 
specific motions into account and the solution will contain only the differences 
t?i - co, i >2 - w and, consequently, one of the angles ^ 1( ^ 2 , co will be arbitrary. 

Thereafter, from the condition (29) for the coincidence of the time spent 
in motion prior to reaching the final point, it is easy to find, for example, the 
angle ^ 


i>i = [ ( o; - 0 )p?p! + coiP 2 - co 2 p? + fpi](p| - Pi) 1 , (42) 

where ip/K is the time spent in motion along the optimum orbit; f = t? 2 - i?i is the 
difference between the true anomalies of the starting and finishing points. 

Let us consider a problem about a flight without taking specific motions 
into account. It follows from the equations ( 36) , ( 37) that 

\i sin (i?i - co) = 0 , \ 2 sin (i? 2 - co) = 0 . (43) 

The last equations can be satisified by using three methods, since the conditions 
\i = 0, \ 2 = 0 contradict the equation (39). 

Method I. Assume that 

sin (i?i - co) = 0 , sin (i >2 - co) = 0 . (44) 

The only solution in this case (with an accuracy of up to an arbitrary choice of 
the angle ^ of the starting point) is a Homan ellipse in geometrical elements. 

For it, we have 


11 



(45) 


4 


P? + P2 


Pi “ P2 


= ± 


4T( P ? 


. 2 
+ P2 


Xi = 


pWi - 5 - PiP 
4p 4 


„ _ Pi(pWl - 6 - p) 

\ 2 — 

4p 4 


AU = ^ I pi n/ 1 - 5"- pi 

P 


(46) 

(47) 


The upper symbols are for flights to orbits having a large radius 

r 2 > r i > P 2 < P < Pi , t?i = co , = co + 7r , 4> = 0 at p < Pl ^1 - 5 , 

$ = 7r at p > p^l - 6 ; (48) 

the lower symbols are for flights to orbits of smaller radius 

r 2 < rj , P 2 > P > Pi , i?i = co - 7r , t? 2 = co , $ = 7 r . (49) 

Method II. We will satisfy the equations (43) , as follows 

sin ( $\ — co) = 0 , A 2 ~ = 0 . (50) 

From the remaining equations of the system of required conditions, we find 

p= PlN fi— T, q = (51) 

2 2 

— co , cos u 2 - co) = P2 ~ Pl(1 ~ , (52) 

P?5 

AU = 0 . (53) 


If 6 < 0.5, a solution exists only at 


1-25 


> r 2 > rj , 41 - 26 Pi < p 2 < Pi , 


(54) 


however, if 6 0.5, the solution is at r 2 > r t . 

The take-off takes place at the pericenter of the flight's orbit, no addi- 
tional consumption of fuel is required and, therefore, these orbits offer a better 
advantage in terms of energy than a Homan ellipse. For flying shell-probes to 
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. reach the obtained orbits, it is sufficient to inflate the shell. Note that at 

r 2 = rj/1 - 26, the obtained orbit will coincide with a Homan ellipse; it will be 
a parabola at 5 = 0.5 and an hyperbola at 6 > 0.5, at which the Sun is located 
in the internal focus. 

Method III . In this case 

Xj = 0 , sin (i? 2 - w) =0 . (55) 


(56) 


(57) 

j Such orbits are possible only at 

i 

o- 2 r! 4 : r 2 < rj , p 2 > p t > op 2 , ( 58) 

| where cr is the only positive root of the equation 

CT 8 + o 6 = 2(1 - 8) . (59) 

I 

The finishing point coincides with the pericenter of the flight’s orbit. At 
r 2 = a^Tx, the obtained orbit coincides with a Homan ellipse. The characteristic 
velocity and the tangent of the angle of inclination of the thrust are determined as 
follows: 

I 

(60) 
(61) 

During the same interval (58) , a Homan-type of a flight is possible (45) -(49) . 

Let us compose the difference between the squares of the characteristic velocities 
along a Homan ellipse and along the new orbit, which is equal to 

kV 2 [rf 

Pi 2 + p! l 


'pT+pT 


- n2( 1 


- 5)J 2 > 0 . 


AU = K'jspi - 2pi - + 2 6 (pi - p?)i 

* Do / 


t $ _ p| tPi - P^ 1 ~ j) 1 s i n Ui ~ to) 
Pi(P? - P$) 


The remaining equations of the system of required conditions will yield 


pWl - <5 



Pi 


= P? - p|U - 6) 

P?pS^1 - 5 


Pi - P2<1 - g) 

COS (* - «) = p 2 [p G _ p 6 (1 _ 6)] ’ 


$■ 2 ~ CO . 
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Consequently, the obtained orbits have a better optimum than a Homan ellipse. 

Therefore, the function (34) assumes its least value on orbits described 
in II and III, provided that the conditions of (54) or of (58) are satisfied, 
respectively; otherwise, the least value will be on a Homan ellipse (45) -(49) . 

NOTE 4. The value of the parameter 6 in many problems is small and in 
such cases the approximate magnitude of the value ct is 


<7 = 



43 

686 


6 2 + 


(63) 


5. Flight Between Low-Eccentricity Orbits 

Assuming that q l9 q 2 are small, let us introduce a small parameter e as 
per the expressions 


qi = q{e » q2 = q&e . (64) 

The equations (27) -(29) and (36) -(40) will be solved in form of series of powers 
of e; we will use a vinculum for the coefficients of the sought series with the 
first powers of e. Similar expansions were obtained by Kirpichnikov 1 for 6=0. 

The solution will be carried out with the assumption that r t , r 2 are not 
located in the zones of (54) or (58) ; i. e. , a Homan ellipse will be taken as 
zero approximation. 

Let us first not take into account the specific motions. In zero-order 
approximation for e, we have only one solution — a Homan ellipse (45) -(49) 

(with an accuracy of up to an arbitrarily chosen one of the angles of ^ , i? 2 , w) . 
The angular distance of the pericenter co is determined by the first approximation 
with the aid of equations (36) , (37) , and (40) , as follows: 


tgcu = 


qlpi 

|\ , 2 Pl 2 - 6_1 

5 p 's/I - 6_ 

sin coi - q2P2^2 si n w 2 

q]pi 

R ^ 2 Pi 2 - 6.- 

_ 1 p n /1 - 6_ 

COS COi - q$> 2^2 cos 


(65) 


The following expressions provide an exact solution of the first order: 


1 Ibid. 
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_ QlPl cos (^1 - w l) + < ltf ) 2 GOS (^2 “ w 2) 

P 4^ ’ 

, ^tlPl ( Q T 2 P) cos (^1 “ w l) + q^Mq ± 2 P) cos (^2 ” 03 2 ) 

p = w 


G d 


f _ ~ a * A - fi* B + C 
a% + (3* + y?,<. 


, i>i = 0) T + A , t>2 ~~ + ® > 


_ 1 pR* - ( q ¥ 2p) S* » pR* - (q? 2p) Sy^ 
A 1 = — » a 2 ~ 


AU 


4p 2 “ 4p z 

K 2 \/ 1 - 6 /pi - p 3 n/i - 5 pj T -7 

’ = ATT r 2 p + qNfl - 6p 


(66) 

(67) 

( 68 ) 


q 1 'p 1 (2p - 2pjN/l - <5 - p 3 ) cos (^ - oq) 
pN /l - a 




$' = ± 


pqVl - 5 (^ - co*) q]p sin (t?i - ag) 
P^Pjn/T^T - p) Pi(pWl - 5 - p) 


(69) 

(70) 


where 


A = ± + £U X, - j^JJsin <* - W ,) , ) 

*ipqLpi p vi - 6J 

► 

B = T q 2 P 2 sin (_&z ~ co 2 ) 

pq ’ ' 

/ 2 2 \ 

= qu— — 7 -^- p’ - PiA-i + q') sin (t?i - Wi) - q 2 * 2 P 2 sin (^2 “ 03 2) > 
\ n 2 ' 


(71) 

(72) 


a* = qiPijjp + Y + X i " ^%%] cos ( * “ Wl) ’ 


/S s jc — q 2 P 2 ^-2 COS (t?2 - k> 2 ) » 

y* = ± qjq cos (i?! - o>i) , 

R* = q 1(1 -—Vos ( - coi) 

\ p2 / 

+ ( Vl - a"+^l - 2\j - 2\ 2 V T (Xi - \ 2 )q' , 

' n 3 / 


(73) 


(74) 
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S* = ± qj cos (^ - coi) + (Xi - ^ 2 ) P* “ “ 5q' . (75) 

The solution of a problem of "energy -optimum" one-thrust flight without 
taking into account the pressure of light 1 with an accuracy of up to the terms of 
the first order of e is obtained with the aid of the expressions (45) -(49) , 

(65) -(75) , provided that it is assumed that 5 = 0. In case the parameter 5 of 
the same order as the eccentricities of the boundary orbits, is 


5 = eS' , (76) 

for this solution, the corrections of the first order which are made necessary 
by the presence of light-pressure, will be 


p' = 0 , q T = 0 
CO* = = 1?2 

2qi l q;>pfpi>p 3 (p 1 - p) 2 sin ( qjp - coi) cos 2 co 

{qiPi [ 8p 3 ( Pi -p ) + P 2 - Pip] cos W! - q^PiP 2 (Pi - p) cos co 2 }* 

x; = -^ 5 ' • = ■ 


AU 1 = + 6 ’ , 4>’ = 0 . 

2p 


1(77) 


Consequently, for a case where the parameter 6 and the eccentricities 
of the boundary orbits are small, the effect of the pressure of light which was 
taken into account was reduced to a turn of the orbit of the flight by a value of 
the first order with respect to an orbit corresponding to an optimum problem in 
which the pressure of light is not taken into account. The size and the shape of 
the orbit and the true anomalies of its starting and final points change only by 
values of the second order. 


Let us proceed with solving the problem with taking into account the 
specific motions. With a zero-order approximation, the only solution is a 
Homan ellipse (45) -(49) which satisfies the equality (41) , while for the function 
ip, we have 


1 Ibid. 


4 > = 


7T 


(p 2 - q 2 ) 8/2 1 - 6 


(78) 
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The starting angle ^ is found from the equality (43) , as follows: 


f + p?^ 2 _i 


t?i = 7T- 


Pl + P2 
. 2 


V 2 -i~ 

7 ^ 


Pl “ P2 


d ~ Pl x co 2 pj - q»iP2 - qpipj 
Pi - P2 


In the approximation of the first order, from the system of required 
conditions we find for p', q' the same expressions (66) , while for t? 2 , and tx>, 
we have: 


, _ A(pfp 2 ^l - d - ppf) - B( pip 2 ^l - <5 - ppf) - Cpfpf\/1 - 6 
” 


PiP 2 (Pi - P 2 ) ~ d 


, _ A^pfp^l -6 - ppf) - B( pfp t \ll - 5 - ppf) - Cpfpf^l - 5 
v2 ~ — — — 

PiP 2 (Pi - P 2 ) - <5 


_ A(p|p 1 N/l - 6 - ppf) - B(pfp 2 s/1 - 6 - ppf) - Cpfpf\/1 - 6 , 
co — — — 

PiP 2 (Pi - P 2 ) 'll - 6 


where 


A _ A _L PWl - d - P , ! -rT _ r, - P - P - O . . 

A \ lP fpqVl - * X 3 . B = B + 7f X 3 > 


p - pWl - d , , 
\ 2 P2pq^l - d 3 


~ 37r(pp' - qq') , 2qJ sin (^ - u>i) 2qf sin (t? 2 - co 2 ) 

v_^ — ■ " ' """" " f ^ 1 "" “ - 


(p 2 - q 2 ) 


5/2 sfT 


the function A and B are computed with the aid of expressions ( 71) , and 


qiPi(— + M - = 8 f j sin (A - wi) + q 2 P 2 * 2 sin (t? 2 - <^ 2 ) 

^3 = ‘ ' 

„ “3 „ “3 

P 2 - Pl 
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6. Condition for Energy-Saving Use of Shell-Probes 
Carrying Useful Loads 

Let us consider in which case it is more profitable (by using less fuel) 
to deliver to a final orbit a useful load with a mass m 0 : with the aid 
of a hollow shell-probe, or by delivering this load directly to a corresponding 
orbit of the flight. In the latter case, the pressure of light can be^disregarded. 
We will assume that both cases require a characteristic velocity U to overcome 
the attraction of the body in space from which, or from the artificial satellite 
of which the start takes place. We will confine ourselves with a case of circular 
boundary orbits, where r 2 > r*. At r t > r 2 , it is more profitable to launch a 
load so, that the resultant of the forces of repulsion by light should be as small 
as possible. 


In a flight in which the pressure of light is not taken into account, a 
Homan ellipse requires a minimum of fuel consumption. The characteristic 
velocity AU^ of the initial thrust is defined as follows: 


AU - K Pi(Pi ~ P) 
r p 


(86) 


The minimum value of the characteristic velocity of a flight which takes into 
account the pressure of light (6 < 0.5), is equal to: 


AU = 0 at 


1-26 


> r 2 > r 4 


AU = £EdjWHl^ P> at r, > - r ' . 

p 1 - 2o 

Let cj be the velocity of the escaping gases; then, as we know, 


(87) 

( 88 ) 


m H „ m - m, + Am 

AU + U = c t In , AU + U = c t In , ( 89) 

nig + m 1 F m 0 

where m^ is the initial mass of the space vehicle at the launching of the useful 

load with the shell-probe; Am is the change in the mass of the fuel during a 
direct launching of the useful load to the orbit; m t is the mass of the shell -probe, 
for which the following is valid: 


nij = 47rR ? yh 


(90) 
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The expressions ( 92) andj 93) make it possible to find the required relationship 
between y, h, x. Pi, P 2 , U. 
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